High-capacity hydrogen storage by metallized graphene 
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First-principles plane wave calculations predict that Li can be adsorbed on graphene forming a 
uniform and stable coverage on both sides. A significant part of the electronic charge of the Li- 
2s orbital is donated to graphene and is accommodated by its distorted 7r*-bands. As a result, 
semimetallic graphene and semiconducting graphene ribbons change into good metals. It is even 
more remarkable that Li covered graphene can serve as a high-capacity hydrogen storage medium 
with each adsorbed Li absorbing up to four Lb molecules amounting to a gravimetric density of 12.8 
wt %. 



Developing safe and efficient hydrogen storage is 
essential for hydrogen economy^ Recently, much ef- 
fort has been devoted to engineer carbon based 
nanostructures 2 ^^ which can absorb H 2 molecules with 
high storage capacity, but can release them easily in the 
course of consumption in fuel cells. Insufficient storage 
capacity, slow kinetics, poor reversibility and high dehy- 
drogenation temperatures have been the main difficulties 
towards acceptable media for hydrogen storage. 

Recently, graphene, a single atomic plane of graphite, 
has been produced^ showing unusual electronic and mag- 
netic properties. In this letter, we predict that metallized 
graphene can be a potential high-capacity hydrogen stor- 
age medium. The process is achieved in two steps: Ini- 
tially, graphene is metallized through charge donation by 
adsorbed Li atoms to its 7r*- bands. Subsequently, each 
positively charged Li ion can absorb up to four H2 by 
polarizing these molecules. At the end, the storage ca- 
pacity up to the gravimetric density of <?d=12.8 wt % is 
attained. These results are important not only because 
graphene is found to be a high capacity hydrogen storage 
medium, but also because of its metallization through Li 
coverage is predicted. 

Our results have been obtained by performing 
first-principles plane wave calculations using ultra-soft 
pseudopotentials»£ We used Local Density Approxima- 
tion (LDA), since the van der Waals contribution to 
the Li-graphene interaction has been shown^ to be bet- 
ter accounted by LDA. Numerical results have been ob- 
tained by using VASP, 9 which were confirmed by using 
the PWSCF codeJ^ A plane-wave basis set with kinetic 
energy cutoff 7i 2 |k + G| 2 /2m = 380 eV has been used. 
In the self-consistent potential and total energy calcula- 
tions the Brillouin zone has been sampled by (19x19x1) 
and (9x9x1) special mesh points in k-space for (2x2) and 
(4x4) graphene cells, respectively. Atomic positions in 
all structures are optimized using the conjugate gradient 
method. Convergence is achieved when the difference of 
the total energies of last two consecutive steps is less than 
1CP 6 eV and the maximum force allowed on each atom 
is less than 10~ 2 eV/A. All configurations studied in this 
work have also been calculated by using spin-polarized 
LDA, which were resulted in non-magnetic ground state. 

Adsorption of a single (isolated) Li atom on the hollow 



site of graphene (i.e. Hl-site above the center of hexagon) 
is modelled by using (4x4) cell of graphene with 1.70 
A minimum Li-graphene distance and with a minimum 
Li-Li distance of 9.77 A, resulting in a binding energy 
of El=1.93 eV. Upon adsorption Li atom donates part 
of the charge of its 2 s state to the more electronegative 
carbon atoms at its proximity. Despite the ambiguities 
in determining the atomic charge, Lowdin analysis esti- 
mates that Li becomes positively charged by donating 
q ^0.35 electrons (but q ~ 0.9 electrons according to 
Bader analysis^). The energy barrier to the diffusion 
of a single Li atom on the graphene sheet through top 
(on-top of carbon atoms) and bridge (above the carbon- 
carbon bond) sites are calculated to be AQ=0.35 eV and 
0.14 eV, respectively. 



Lithium atoms can form a denser coverage on the 
graphene with a smaller Li-Li distance of 4.92 A form- 
ing the (2x2) pattern. Owing to the repulsive interaction 
between positively charged Li atoms, the binding energy 
of Li atom is smaller than that of the (4x4) cell. For 
HI adsorption site (see FigQJa)), the binding energy is 
calculated to be El— 0.86 eV. The binding energies are 
relatively smaller at the bridge and top sites, and are 
0.58 and 0.56 eV, respectively. The binding energy of the 
second Li for the double sided adsorption with H1+H2 
and H1+H3 configurations described in FigJT] (a), are 
El =0.82 and 0.84 eV, respectively. The same binding 
energies for H1+H2, and H1+H3 geometries on the (4x4) 
cell are relatively larger due to reduced repulsive Li-Li in- 
teraction, namely El=1.40 eV and 1.67 eV, respectively. 
The coverage of Li on the (2x2) cell is 0=12.5 % (i.e. 
one Li for every 8 carbon atoms) for HI geometry and 
0=25 % for either H1+H2 or H1+H3 geometries. Metal- 
lic charge accumulated between Li and graphene weakens 
the interaction between Li atoms which are adsorbed at 
different sites of graphene. Further increasing one-sided 
coverage of Li to 0=25% with HI geometry (or two-sided 
coverage to 50 % with H1+H2 or H1+H3 geometries) 
appears to be impossible due to strong Coulomb repul- 
sion between adsorbed Li ions and results in a negative 
binding energy (El ~ - 2.5 eV). On the other hand, the 
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FIG. 1: (Color online) (a) Various adsorption sites HI, H2 and H3 on the (2x2) cell (top panel) and energy band structure of 
bare graphene folded to the (2x2) cell (bottom panel), (b) Charge accumulation, Ap + , calculated for one Li atom adsorbed 
to a single side specified as HI (top) and corresponding band structure, (c) Same as (b) for one Li atom adsorbed to Hl-site, 
second Li adsorbed to H3-site of the (2x2) cell of graphene. Zero of band energy is set to the Fermi energy, Ef- 



total binding energy of all Li atoms adsorbed on a(2x2) 
cell with the H1+H2 (H1+H3) geometry corresponding 
to 6=25 % is 3.23 eV (3.12 eV) higher than that of Li 
atom adsorbed on the (4x4) cell with the same geometry 
corresponding to 8=6.25 %. Hence, since the cluster for- 
mation is hindered by the repulsive interaction between 
the adsorbed ions, a stable and uniform Li coverage on 
both sides of graphene up to 9 = 25 % can be attained. 



The charge accumulation and band structure calcu- 
lated for the HI and H1+H3 adsorption geometries are 
presented in FigJT] (b) and (c), respectively. Isosurfacc 
plots of charge accumulation obtained by subtracting 
charge densities of Li and bare graphene from that of 
Li which is adsorbed to graphene, Ap + , display positive 
values. As a result of Li adsorption, the charge donated 
by Li is accumulated between Li and graphene and is 
accommodated by 2p7r*-bonds of carbons. The empty 
7r*- bands become occupied and eventually get distorted. 
Occupation of distorted graphene- 7r* bands gives rise to 
the metallization of semimetallic graphene sheets. By 
controlled Li coverage one can monitor the position of 
Fermi energy in the linear region of bands crossing at the 
K-point of the Brillouin zone. Metallization is also im- 
portant for zigzag and armchair graphene nanoribbons, 
since both are semiconductors with their energy gaps de- 
pending strongly on the widths of these ribbons. 12 Seg- 
ments of these ribbons metallized by Li adsorption may 
be interesting for their electronic and spintronic applica- 
tions. For example, a junction of two nanoribbons with 



and without Li adsorbed segments can serve as a Schot- 
tky barrier. 

Sodium, a heavier alkali metal, can be bound to 
graphene with £^=1.09 eV at HI site. However, the 
energy difference between top, bridge and HI sites are 
minute due to relatively larger radius of Na. Upon ad- 
sorption, graphene and graphene nanoribbon are met- 
allized. Nevertheless, Na is not suitable for hydrogen 
storage because of its heavier mass and very weak bind- 
ing to H2 molecules. Two dimensional BN-honeycomb 
structure, being as a possible alternative to graphene, 
has very weak binding to Li (~ 0.13 eV) and hence it is 
not suitable for hydrogen storage. 

The absorption of H2 molecules by Li + graphene in 
HI, H1+H2 and H1+H3 geometries. A summary of our 
results about the H 2 absorption are presented in Fig. 
[21 The binding energy of the first absorbed H2, which 
prefers to be parallel to graphene, is generally small. 
However, when two or more H2 molecules are absorbed 
by the same Li atom, the binding geometry and mech- 
anism change and result in a relatively higher binding 
energy. All H2 molecules are tilted so that one of two H 
atoms of each absorbed H2 molecules becomes relatively 
closer to the Li atom. A weak ionic bond forms through 
a small amount change (> 0.1 electrons) transferred from 
Li and graphene to nearest H atoms of absorbed H2 
molecules. At the end, H atoms receiving charge from 
Li becomes negatively charged and the covalent H2 bond 
becomes polarized. Weak ionic bond, attractive Coulomb 
interaction between positively charged Li and negatively 
charged H and weak van der Waals interaction are re- 
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FIG. 2: (Color online) Adsorption sites and energetics of Li adsorbed to the (2x2) cell of graphene and absorption of H2 
molecules by Li atoms. El is the binding energy of Li atom adsorbed to Hl-site, which is a minimum energy site. For H1+H2 
or H1+H3 configuration corresponding to double sided adsorption, El is the binding energy of second Li atom and El is the 
average binding energy. For HI, H1+H2 and H1+H3 configurations, Ei is the binding energy of the first H2 absorbed by each 
Li atom; E n (n=2-4) is the binding energy of the last n th H2 molecule absorbed by each Li atom; E n is the average binding 
energy of n H2 molecules absorbed by a Li atom. Shaded panel indicates the most favorable H2 absorbtion configuration. 



sponsible for the formation of mixed weak bonding be- 
tween H 2 molecules and Li+graphene complex. Here the 
bonding interaction is different from the Dewar-Kubas 
interactioniS found in H2-Ti+Cgo or carbon nanotube 
complexes^ As the number of absorbed H2, n, increases, 
the positive charge on Li as well as the minimum distance 
between H2 and Li slightly increases. No matter what the 
initial geometry of absorbed H2 molecules would be, they 
are relaxed to the same final geometry presented in Fig[2] 
for any given n. We found no energy barrier for a H2 
molecule approaching the absorbed H2 when n < 4. Note 
that the dissociative absorption of H2 molecules do not 
occur in the present system. The energy barrier for the 
dissociation of H2 near Li to form Li-H bond is ~ 2 eVJ^ 
Moreover, dissociation of H2 to form two C-H bonds at 
the graphene surface is energetically unfavorable by 0.7 
eV. 

Maximum number of absorbed H2 per Li atom is four, 
and the maximum gravimetric density corresponding to 
H1+H2 geometry at 0=25 % coverage is 5^=12.8 wt 
%. This is much higher than the limit (gd=(> wt %) set 
for the feasible H2 storage capacity. Note that only for 
n = 4, H1+H2 absorption geometry has slightly lower 
energy than H1+H3 geometry. This is a remarkable re- 
sult indicating another application of graphene as a high 
capacity storage medium. Here Li+graphene complex 
is superior to Ti+Cgo or carbon nanotube complexes 
since Li is lighter. Even though graphene by itself is 
stabl o 15 i 16 , more stable form is obtained by Li adsorption 



on graphene due to strong Coulomb repulsion between 
adsorbed Li atoms. Moreover, Li covered graphene is re- 
sistant to clustering of adsorbed Li atoms. Earlier, Dur- 
gun et al£ has predicted that ethylen+Ti complex can 
store H2 up to <7d=l4.4 wt % per molecule. Later, their 
results have been confirmed experimentally. 17 We believe 
that hydrogen storage by the Li covered graphene is in- 
teresting, since it may not require encapsuling and hence 
may yield even higher effective ga- 
in conclusion, two crucial features of Li covered 
graphene revealed in this paper may be of technological 
interest. These are high metallicity and high hydrogen 
storage capacity of graphene functionalizcd by Li atoms. 
Graphene nanoribbons metallized through adsorbed Li 
atoms can be used as interconnects between graphene 
based spintronic devices. Graphene functionalized by Li 
can also serve as a medium of hydrogen storage. As far 
as efficiency in storage is concerned, graphene may be 
superior to carbon nanotubes, because its both sides are 
readily utilized. Cell configurations formed by different 
junctions of graphene^ functionalized by Li atoms are 
expected to yield higher surface/ volume ratio and hence 
to provide efficient H2 storage in real applications. 
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